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Abstract 

The free energy of U(N) gauge theory is expanded about a center-symmetric 
topological background configuration with vanishing action and vanishing 
Polyakov loops. We construct this background for SU(N) lattice gauge the- 
ory and show that it uniquely describes center-symmetric minimal action 
orbits in the limit of infinite lattice volume. The leading contribution to the 
free energy in the 1/N expansion about this background is of O(N ) rather 
than 0(N 2 ) as one finds when the center symmetry is spontaneously broken. 
The contribution of planar 't Hooft diagrams to the free energy is 0(1/N 2 ) 
and sub-leading in this case. The change in behavior of the diagrammatic 
expansion is traced to Linde's observation that the usual perturbation se- 
ries of non-Abelian gauge theories suffers from severe infrared divergences j^. 
This infrared problem does not arise in a center-symmetric expansion. The 't 
Hooft coupling A = g 2 N is found to decrease cx l/ln(iV) for large N. There 
is evidence of a vector-ghost in the planar truncation of the model. 
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1 Introduction 



Confinement can be denned as the absence of asymptotic states in non-trivial 
multiplets of the global gauge group. Since the number of singlet states does 
not increase proportional to A, the free energy of a U(N) gauge theory 
in a confining phase should be of order iV°[2j. Perturbatively, the adjoint 
multiplet of gauge bosons and the fundamental fermion multiplets contribute 
to the free energy density of U(N) gauge theory in 0(N 2 ) and O(N). A 
direct application of 't Hooft's 1/N expansion [3J H] apparently also gives a 
free energy density of order A 2 even at low temperatures. Using the Dyson- 
Schwinger equations of the lattice and Migdal's factorization condition for 
planar diagrams, Gocksch and Neri[5 j on the other hand found that the free 
energy density in the confining phase does not depend on the temperature 
at A = oo. 

A leading contribution to the free energy of U(N) gauge theory of order 
A 2 is not compatible with the result of Gocksch and Neri[5|. It is more 
reasonable to assume that the coefficient of the A 2 -term in the 1/A-expansion 
of the free energy vanishes in the confining phase and that the model defined 
by planar diagrams is a topological theory without dynamical degrees of 
freedom. It turns out that this is the case in the A —>■ oo limit only. The 
factorization for large A in the confining phase implies that U(N = oo) is 
described by a matrix model that depends on space-time parametrically[UJ 
i only. 

The confining phase of pure Yang-Mills models is characterized by a global 
center symmetry. This symmetry also is essential in the formulation of re- 
duced models[7l EJ E| at N — oo. The objective here is to construct an 
1/N expansion that preserves this center symmetry in every order. Since the 
pure Yang- Mills action is invariant, one can achieve this by expanding about 
a center-symmetric topological field configuration. 

We show the absence of contributions to the free energy proportional 
to A 2 and A when SU(N) gauge theory is expanded about such a center- 
symmetric orbit. In this expansion, planar 't Hooft diagrams contribute to 
the free energy in order 1/N 2 . The present analysis systematizes and extends 
the result of Gocksch and NerijS] in several ways. The center-symmetric 
1/N expansion is possible for all A. It not only gives the order of planar 
contributions to the free energy but also of higher genus 't Hooft diagrams. 
Although fields in the fundamental representation explicitly break the center 
symmetry, there are no contributions to the free energy of order A in this 
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expansion - planar diagrams with a single fundamental color loop contribute 
in 0(1/N 3 ). The leading temperature dependent contributions to the free 
energy are of order N°. These non-planar contributions survive the large N 
limit - as might be expected if the masses of asymptotic singlet states have 
a finite limit 4.. 

In the topological sector with vanishing instanton number that interests 
us here, the classical action of a gauge theory vanishes at field configurations 
with minimal action. However, this classical field is not necessarily a pure 
gauge configuration. Since the local curvature of the configuration vanishes, 
the possibly nontrivial gauge invariant quantities are non-contractible Wilson 
loops. [These non-contractible loops in general are sensitive to global sym- 
metries of the action and thus can distinguish different phases of the model.] 
At a finite temperature T and infinite volume V, configurations with van- 
ishing curvature are characterized by their Polyakov loops, non-contractible 
Wilson loops in the Euclidean temporal direction. 

Specifically, consider the Polyakov loop of an SU(N) gauge theory at 
finite temperature T with periodic boundary conditions for the connection, 



In Eq. (JTJ V denotes ordering of the exponential along the path and is the 
gauge connection in the fundamental representation. On the lattice, U(x) is 
the ordered product of the links in the periodic temporal direction, beginning 
with the link at x. One can choose a gauge in which V^x, r) does not depend 
on the Euclidean time r and is diagonal. On the lattice this may be achieved 
in three steps: one first uses the gauge freedom to set all temporal links 
apart from those on the £4 = time slice to unity. [The nontrivial temporal 
links of this representative configuration then are the Z7(x, 24 = 0) of the 
Polyakov loop.] One next uses time-independent gauge transformations to 
diagonalize the remaining nontrivial temporal links. Since the permutation 
group is a subgroup of SU(N), the phases in addition can be ordered so that 
the temporal links of the x^ = time slice are of the form, 



£(x) = TrU(x) with U(x) = Pexp % 




N 



U(x,x 4 = 0) 



= diag(e iei(x) , . . . , e^ (x) ), with e i( x ) = > 



and - 7T < 0i(x) . . . < 0j(x) < j+ i(x) . . . < ^(x) < vr ,(2) 
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and all other temporal links are unity. The Abelian invariant subgroup of 
the configuration is enhanced to a non- Abelian one when some of the phases 
in Eq.Q are degenerate. The corresponding continuum configuration in 
this case may have a non-trivial monopole number 10]. [Since all lattice 
configurations are contractible, the usual topological classification of smooth 
continuum configurations cannot be used, but degenerate configurations that 
are invariant under a non- Abelian subgroup of SU (N) can also be found on 
the lattice.] 

One finally may use time-dependent Abelian gauge transformations to 
evenly distribute the U(x, 0) of Eq.(j2J) in temporal direction. In the contin- 
uum limit, the resulting configuration corresponds to a temporal component 
of the connection V^x) that does not depend on the Euclidean time x 4 and 
is Abelian. 

The perturbation series can be constructed about any configuration of 
minimal classical action. Although such a configuration generally will not 
correspond to a minimum of the effective action, the perturbation series 
nevertheless yields some information about the configuration space in its 
vicinity. We will see that the confining phase extends to arbitrary small 
values of the 't Hooft coupling for sufficiently large N. That the effective 
coupling may become weak in the confining phase at sufficiently large N was 
previously observed by exploiting the analogy with string theory. We 
show that the perturbative analysis leads to the same conclusion. 

In the topologically trivial sector, the local curvature of a minimal action 
orbit vanishes. The previous construction implies that the temporal links 
of periodic lattice configurations with minimal Wilson action can be chosen 
Abelian and constant across the whole lattice. Since every plaquette-action 
of a minimal action configuration vanishes and all temporal links apart from 
those on a particular time slice can be set to unity by a gauge transformation, 
the spatial links of a minimal action configuration do not depend on time 
in such a gauge. Periodicity of the configuration in time then requires that 
the eigenphases of two spatially adjacent temporal Abelian links are the 
same: since all plaquette-actions vanish we must also have that ga = a'g, or 
gag' = a' for two equal spatial links g G SU(N) and two adjacent temporal 
links a and a' on the £4 = time-slice. The previous procedure shows that a 
and a' can be chosen to lie in the Abelian subgroup of SU(N). a and a' thus 
are the same up to a permutation of their eigenphases. Taking into account 
that the eigenphases have been ordered, one concludes that a = a' in this 
particular gauge. 
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All temporal links on the X4 = time slice of this representative of an 
orbit with minimal Wilson action thus are Abelian and the same - all other 
temporal links are unity. We in particular have that minimal action config- 
urations of a time-periodic S{7(iV)-lattice are characterized by a Polyakov 
loop that does not depend on the chosen spatial point. A spatially constant 
Abelian gauge transformation can be used to evenly distribute the temporal 
links of the X4 = time-slice in temporal direction. One thus obtains a rep- 
resentative of any orbit with minimal Wilson action that is described by a 
temporally and spatially constant Abelian connection V4. When none of the 
eigenphases of the temporal links are degenerate, spatial links in this gauge 
also have to be in the Abelian subgroup and do not depend on Euclidean 
time. 

2 Topological configurations and center sym- 
metry 

The minimal action configurations of SU(N) are further characterized by 
their transformation under a global Zn symmetry of the Wilson action. This 
so-called center symmetry is generated by multiplying every temporal link 
on a particular time slice by an element of the center of SU(N) - possi- 
bly followed by a (periodic) gauge transformation of the configuration. This 
transformation multiplies the Polyakov loops of any configuration by a root 
of unity, but does not change the Wilson action. The center symmetry there- 
fore maps minimal (Wilson) action configurations onto themselves. It allows 
to distinguish between minimal action orbits that are invariant under this 
discrete global symmetry and those that are not. 

2.1 The center-symmetric topological configuration 

Since any Polyakov loop is multiplied by a root of unity, an orbit is center- 
symmetric only if its Polyakov loops vanish. The N eigenphases of U{x) G 
SU (N) in Eq.([T|) therefore sum to zero and their product is detU(x) = 1. The 
discussion in the introduction shows that one may choose Z7(x, 0) constant 
and in the Abelian subgroup. The constant 6^ on the X4 = time slice of 
such a center-symmetric minimal action configuration thus are, 

Of (x) = tt(2j - N - I) /N , for j = 1, 2, . . . , N . (3) 
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A center transformation simply permutes the phases in Eq.fjHJ) and the 
previous ordering can be restored by a time-independent SU (N) gauge trans- 
formation (of which the permutations are a subgroup). The fact that the 
eigenphases in Eq.Q are equidistant was recently exploited to define an 
order parameter for the center-symmetric phase [T^]. 

None of the eigenphases of a center-symmetric configuration with minimal 
action are degenerate. The spatial links therefore do not depend on time and 
are Abelian as well. On a lattice that is periodic in every direction they 
can in fact be chosen Abelian and constant. To see this, one may proceed 
as follows. Using time-independent Abelian gauge transformations only, all 
(already time independent) spatial links in ^-direction apart from those on 
the 2 3 = slice may be set to unity. This time-independent Abelian gauge 
transformation does not change the temporal links. Since this is an Abelian 
minimal action configuration on a lattice that is periodic in 23, the links in 
the 23-direction on the 23 = slice in fact must all be equal. The remaining 
Abelian links in the 22- and 21-directions at this stage do not depend on 
23 (nor on 24). Using an Abelian gauge transformation that depends on 
23 only, the links in 23-direction on the 2 3 = slice can be distributed 
evenly in the 23-direction. The result is a gauge equivalent configuration 
with constant Abelian links in 24- and 23- directions and Abelian links in 
2 2 - and x\- directions that do not depend on 2 3 nor on 24. The procedure is 
repeated with 24 and 23-independent Abelian gauge transformation to also 
make the links in 22-direction constant (links in 21-direction at this point do 
not depend on 22,23 nor 24). Abelian gauge transformations that depend 
only on x\ can finally be used to obtain a configuration with Abelian links 
in each direction that do not depend on space or time. 

In general there are inequivalent center-symmetric minimal action orbits 
that differ in the eigenphases of the spatial links. However, this distinction is 
critical at finite volume only. The above construction implies that the phases 
of the constant spatial links of the final configuration can be chosen to all fall 
in the interval (—n/L, tt/L], where L is the spatial lattice dimension in lattice 
units. In the limit L — > 00, the spatial links of the configuration all tend 
to unity. The arbitrarily small deviations from unity can only be observed 
by non-contractible Wilson loops that wrap around the whole spatial extent 
of the lattice. These are not observables in the infinite volume limit and a 
center-symmetric orbit of minimal action in this sense is unique. 

In the infinite volume limit at a finite temperature T any center-symmetric 
orbit with vanishing curvature can be represented by a constant Abelian con- 
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nection. Using Eq.fjSJ) and the previous observation that spatial links of this 
representative tend to unity for large spatial volume, this center-symmetric 
Abelian background connection is, 

gV} 0) =ga 4 = Tdiag(^ 0) , . . . , < } ) 



2nT diag 



-N +1 -N + 3 JV-3 JV- 1 



2N 2N 1 2N 1 2N 

gV^ =g ai = 0, for % = 1,2, 3. (4) 

Quadratic fluctuations about the center-symmetric configuration of tem- 
poral links have been considered previously (see for instance O EH] and (for 
A" — ► oo) in [IT]). We here will examine the perturbation series about the 
configuration Eq.pjl to all orders in the 1/N expansion of the free energy. 

2.2 Topological configurations that break the center 
symmetry 

If any Polyakov loop of a topological configuration does not vanish, it nec- 
essarily belongs to a multiplet of minimal action configurations. If A" is not 
prime, the configuration may break a subgroup of Z^ only. However, the 
flat connection V4 = breaks the Z^-group completely. It is one of the N 
Abelian configurations of the form, 

# 4 (9) = ^diag((l-A0,l,...,l,l); q=l,2,...,N (5) 

These configurations have degenerate eigenphases and one cannot argue that 
the spatial links of such a minimal action configuration are Abelian. Con- 
trary to the center-symmetric case, it is not clear that the index q uniquely 
identifies a minimal action orbit in the infinite volume limit. 

The configurations of Eq.(jSJ) have been studied extensively |IH]. They 
correspond to minima of the free energy at high temperatures T when cor- 
rections proportional to the coupling g 2 {T) are negligible. [However, the 
homogeneous vacua of Eq. (j3J) do not solve the infrared problem of the high- 
temperature expansion observed by Linde[l a - the high-temperature phase 
probably [in] can be described by domains of such vacua with different index 

Minimal action configurations that break the center-symmetry to a sub- 
group of Zn also can be constructed for non-prime N. They could play a 
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role in the (perhaps rather complex) phase structure of an SU(N)-mode\ 
with non-prime N. Minimal action solutions that break the Z(iV)-symmetry 
correspond to perturbative minima of the free energy. They are not the 
minima of the free energy in a center-symmetric (confining) phase. 

3 Large N expansion in a center- symmetric 
background 

We are interested in the expansion of the free energy of a [/(iV)-model at 
finite temperature for large values of N. We shall argue that the model is in a 
confining phase as long as the center-symmetric background is stable. [More 
specifically, the free energy density F of a U(N) gauge theory expanded 
about the center-symmetric background is O(N ) rather than (D(N 2 ) and 
0(N) as one expects when asymptotic states form multiplets of the adjoint, 
respectively fundamental, representation of SU(N).] 

The center-symmetric background of Eq. (J3J) is a maximum of the 1-loop 
free energy[T3]. whose minima are at the configurations of Eq.([5|) that spon- 
taneously break the center symmetry. It was recently found that the 
non-perturbative contribution from calorons|14j. can make the minima of 
the 1-loop free energy unstable at low temperatures. Near the deconfine- 
ment transition, calorons with non-trivial holonomy have been observed by 
cooling SU(2) and SU(3) lattice configurations |To] . However, classical solu- 
tions of finite action generally are suppressed in the limit of large iV (g 2 N 
finite). A semi-classical mechanism for restoring the center symmetry thus 
appears unlikely at large N. 

Lattice studies at relatively small N see a distribution of values for the 
Polyakov loop in the confining phase, rather than a strong concentration near 
£(x) = 0. This can also be seen by studying the strong coupling expansion in 
a gauge where all temporal links except those on the £4 = time slice are set 
to unity. To leading order, the measure for the eigenphases 0j(x) of Eq.© 
in this case is given by the Vandermonde determinant of the eigenvalues of 
the non-trivial temporal link, 



This measure is gauge invariant and vanishes when any two eigenphases coin- 




i=l 



i>j 
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cide 2 . It is maximal when the iV phases are evenly distributed over the circle 
[0, 2tt]. For small values of N, the dependence on the eigenphases of Eq.fJUJ) 
is rather weak. However, for iV ~ oo the support of the measure Eq.© be- 
comes restricted to the immediate vicinity of the configuration Eq.fjSJ). The 
strong coupling limit of U(N ~ oo) lattice gauge theory thus is an example 
for the more general conjecture J7j that fluctuations are suppressed at large 
N. 

Since its expectation vanishes, the usual factorization argument fails for 
the Polyakov loop in a center-symmetric phase. However, at strong coupling 
one can explicitly show[18 j that the distribution of C = Till = exp(i9j) 
converges to a standard normal in the limit N — > oo. The fact that the 
variance of the Polyakov loop does not grow with N implies that the stan- 
dard deviation of the eigenphases is of order 1 / y/~N only. This also is ap- 
parent from Eq.(jB|). < C 2 >= 1 furthermore is consistent with usual 1/N 
counting, which suggests that the standdard deviation of the eigenphases is 
(9(1/ ^/(iV)) in the confining phase even when the strong coupling limit does 
not apply. 

The orbit described by Eq. (jlj) is invariant under spatial translations and 
rotations and minimizes the free energy density of lattice gauge theory in 
the strong coupling limit. It also is the only center-symmetric candidate for 
a perturbative vacuum orbit. Since fluctuations of the temporal links are 
expected to be small at large N, we now consider the perturbative expansion 
about the background of Eq. (j3J) at large N. 

The Euclidean time derivative of a minimally coupled field in a non- 
trivial representation of the group occurs through the covariant derivative 
only. In the background of Eq.pjl the time derivative of a field $ of the 
adjoint representation thus is replaced by, 

d&l - V^l = d^t + ig[a 4 , m = d^l + —(a - b)<S> a b ; a, b = 1, . . . , N . 

(7) 

The time derivative of fields \1/ in the fundamental representation is sim- 
ilarly replaced by, 

d 4 V a -> V A ^ a = d^! a +igaf^ b = d 4 ^ a +^—(a ^)^ a ; a = l,...,N . 

N 2 (8) 

2 Degenerate configurations that are invariant under a non-Abelian subgroup of U(N) 
thus have vanishing weight at strong coupling. 
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Physical correlation functions are colorless. All color indices are summed 
over. At any finite temperature and for any N Eq. ((7j) and Eq. JSJ) imply that 
we can associate a discrete "color momentum", 



*(a) 



2nT 



( a -(jV + l)/2); a = l,..., 



(9) 



iV 



with a color index of the fundamental representation. For sufficiently large 
iV one is tempted to replace sums over color indices by integrals and neglect 
the error due to the fact that £(a) only takes discrete values, 



Note that color momentum is in the compact interval [— ttT, ttT] that does not 
depend on N. Loop integrals over color momentum do not induce new UV- 
divergences. [Something rather similar occurs in solid-state physics where 
momenta are restricted to a single Brillioun-cell - the associated space is 
an infinite (periodic) lattice of points. There are no UV-divergences in this 
case, since the smallest distance is the lattice spacing.] The factors of N 
(see Eq. fllOjl ) from the loop integrals over color momentum can almost all 
be absorbed by redefining the coupling, 



Contrary to the 1/iV-expansion in the broken phase, the reduced coupling A 
does depend on N. We argue in section 4.2 that the remaining dependence 
is logarithmic only. 

In his seminal work on large iV|3] 't Hooft has shown that the contri- 
bution of a connected vacuum diagram in ordinary perturbation theory is 
proportional to a power of N that depends on the difference in the number 
of color- and momentum- loops of the diagram only. Using 't Hooft's double- 
line notation one obtains a topological expansion of U(N) gauge theory in 
terms of the genus of perturbative diagrams. For a background configuration 
of Eq.fjSJ) that breaks the center symmetry, the topological expansion coin- 
cides with an expansion in powers of 1/N (in powers of 1/N 2 when there are 
no fundamental representations). In the broken phase, contributions to the 
free energy of leading order in N are given by planar 't Hooft diagrams that 
have the topology of a two-sphere, 5*2. 




(10) 



g 2 N - A . 
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Some of the characteristics of the usual 1/N expansion are retained by 
an expansion of the model about the center-symmetric but ^-dependent 
configuration of Eq.pjl. Since the background is diagonal in color, one can 
still follow the color flow using 't Hooft's double-line notation (see Fig. 1). 




Fig.l: 'T Hooft's double-line notation for the "dressed" gluon propagator: inser- 
tions of the diagonal Abelian background of Eq. @ (depicted by shaded squares) 
do not change the color of a line ! 

One therefore still has a topological expansion in the genus of the 2- 
dimensional surfaces described by 't-Hooft diagrams. However, this topo- 
logical expansion in general no longer coincides with an expansion in 1/N. 
If one could neglect the error due to the discreteness of color momentum 
in Eq.([10p. each color loop indeed would contribute a factor of iV only. One 
then reaches the same conclusions about the order of a diagram as in the 
broken case. However, due to the discretization error, diagrams of a given 
genus in the topological expansion may also contribute to higher orders of 
the 1/N expansion. The genus of a diagram thus only gives the lowest (su- 
perficial) order in the 1/N expansion to which it may contribute. This has 
some interesting consequences for the 1/N expansion of the free energy. 

We show below that the contribution of planar diagrams to the free energy 
of a U(N) gauge theory without fundamental fields is of order 1/N 2 in the 
center-symmetric background of Eq.pjl. [Although they did not specify the 
order in 1/N, Gocksch and Neri[S] also found that planar diagrams do not 
contribute at N = oo.] The leading gluonic contribution to the free energy 
is of order A^° and given by 't Hooft diagrams with the topology of a torus. 
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3.1 Planar U(N) at finite temperature 

The flow of color momentum in planar diagrams is closely associated with 
that of ordinary momentum. Consider gluonic (vacuum) diagrams without 
external legs in the double-line notation of 't Hooft[3j. The number of mo- 
mentum loops, L p , of a vacuum diagram with E gluon propagators and V 
interaction vertices is, 

L P = E-V + 1 . (12) 

If the diagram is "planar", it has the topology of a 2-sphere 3J, 5*2. Gluon 
propagators are the edges of cells and the Euler number, x, of a diagram is, 

X = V-E + L C . (13) 

Here L c is the number of faces, V is the number of vertices and E is the 
number of edges of the complex. 'T Hooft's double-line notation shows that 
the number of faces, L c , is just the number of independent traces over fun- 
damental color indices, that is the number of loops over color momentum. 
Eq.flnj) and Eo.fTS]) with x(S 2 ) = 2 imply that, 

L c = L p + 1 for planar vacuum diagrams . (14) 

As indicated in Fig. 2a, the loops over ordinary momentum can be chosen 
to coincide with the color traces in planar diagrams. One can enforce ordinary 
momentum conservation at each vertex by writing the momentum of a gluon 
propagator as the difference of two loop momenta associated with each face 
of the oriented cells the propagator is an edge of. In vacuum diagrams one 
ends up with just as many loop momenta as color traces. However, one of 
these loop momenta amounts to an overall translation of all other momenta 
and is redundant. We again arrive at Eq.([14p. 
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Fig. 2: a) A typical planar gluonic vacuum diagram that superficially is of order 
N 2 . The flow of color and of ordinary momentum on each of the faces a, (3, ■ ■ ■ 
is given by composite loop momenta 1^,1^.... The trace over color for the 
perimeter loop results in a factor of N. b) A typical planar vacuum diagram 
with one fundamental loop that superficially is of order N. The flow of color and 
ordinary momentum again is captured by composite loop momenta but there is 
no trace over color only. 

This association between color- and momentum- loops in planar diagrams 
can be exploited. In equilibrium at finite temperature, gluons are periodic 
fields in Euclidean time with period 1/T. Their Matsubara frequency u n 
therefore is an integer multiple of the fundamental frequency 2ttT, 

n = 2nTn, nGZ. (15) 

We may enforce momentum conservation at a vertex (also at fermion-gluon 
vertices) by writing the Matsubara frequency of a gluon as the difference of 
the temporal components of two half-integer loop momenta associated with 
the faces (say a and (3) it is an edge of, 

0n-m = ki a \n)-kf\m) = 27rT((n+l/2)-(m+l/2)) = 2irT(n-m) . (16) 

For a planar vacuum diagram the loop momenta can be chosen to run 
along the color loops. We thus can combine the time-component of loop 
momentum fc^ with the color momentum to the temporal component 
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of a single composite loop momentum l^ 1 , 

l^\n, a) = k[ a \n) + £ (a) (a) = ^-(Nn + a- 1/2); n G Z, a = 1, . . . , N . 

(17) 

The conservation of the time component of ordinary loop momentum and 
of color at a vertex, thus is equivalent to the conservation of the integer 
j = Nn + a, i.e. the time component of composite loop momentum Z 4 . Note 
that the sum over the temporal component of composite momentum extends 
over all half-integers and that the temperature effectively is T/N in planar 
U(N). In purely gluonic planar vacuum diagrams, every summation over a 
composite loop index apart from one (the "peripheral" color loop) is accom- 
panied by a factor of T. For the peripheral loop of a gluonic planar vacuum 
diagram the summation is over color only. It amounts to a translation of all 
other composite loop momenta by a half- integer between 1/2 and N + 1/2. 
This changes all other summations over half-integer composite loop momenta 
to summations over integer composite loop momenta and in addition yields 
an overall factor of N [since the expression for the diagram in fact does not 
depend on finite shifts of all composite loop momenta by integer multiples of 
2ttT/N). 

A planar gluonic vacuum diagram with L p momentum loops is of per- 
turbative order (g 2 ) Lp ~ l and is proportional to a factor NT Lp due to L p 
summations over integer composite loop momenta and the trace over color 
of the peripheral loop. For the background of Eq.(J3]), the regularized 3 planar 
contributions to the free energy density, Fg 2 (T), scale as, 

F S2 (T,g 2 ,N) = NTf S2 (T/N,X) = T 4 /iV 2 / S2 (A) . (18) 

After the UV-regularization is removed /(A) is a dimensionless function of 
the reduced physical coupling A(T/A), where A is the appropriate asymptotic 
scale parameter of the renormalization scheme (see section 4.2). 

Eq. ffTHj) shows that there is no contribution of order A^ 2 to the free en- 
ergy from planar diagrams in the center-symmetric background of Eq.Q. 
This implies the absence of asymptotic states in the adjoint representation 
of the group, that is of (constituent) gluons, in center-symmetric U(N). The 

3 The superficially quartic ultraviolet divergence of the free energy can be reduced to 
the superficially quadratic divergence of the specific heat. It then is sufficient to regulate 
the spatial integrals (see section 4.2). The severe infrared divergences of perturbation 
theory observed by Lindeppj are absent in the present case (see section 6). 
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free energy of the model otherwise would have to be proportional to N 2 , the 
degeneracy of such a multiplet. The result also eliminates the possibility 
of asymptotic states in higher dimensional representations. Eq. (jl8|) sug- 
gests that the leading contribution to the free energy of gluonic and center- 
symmetric U(N) is of orderiV and given by diagrams with the topology of 
a doughnut T 2 . 

Although this is more or less what one would expect for the confining 
phase of the model, some omissions and apparent contradictions have to be 
addressed. Any explicit calculation requires the specification of a gauge and 
an appropriate regularization procedure. We have to show the existence of 
a gauge that is compatible with the background of Eq.(jU) and does not in- 
validate the previous argument. We also still have to verify Eq.(|18|) for the 
(planar) contribution to the free energy of U(N) of order A . This "1-loop" 
contribution to the free energy is a Casimir energy that does not correspond 
to an evaluation of vacuum diagrams like those discussed above. The follow- 
ing sections support the above argument in important ways. 

4 Gauge Fixing and Renormalization 

4.1 Background Gauge 

The background configuration of Eq. (j3J) is in the maximal Abelian subgroup 
of U(N) and a crucial point of the previous argument was that all fields 
couple minimally to it. Covariant Maximal Abelian gauges (MAG) satisfy 
this requirement and furthermore can be defined [21! on the lattice 4 . The 
Abelian Ward Identity of MAG implies that the background ga^ does not 
renormalize in these gauges [21]. It therefore is sensible to set this background 
connection proportional to the physical temperature T in MAG. 

However, the fact that MAG distinguishes between diagonal and off- 
diagonal components of the connection gives rise to additional vertices at 
which the color flow is constrained. In diagrams containing such vertices, 
not all color loops are independent. This leads to apparent modifications of 
the 1/A-expansion and complicates the 1/iV-counting considerably: due to 
cancellations, gauge-invariant combinations of diagrams can be of different 
order in 1/N than the connected diagrams are individually. 

4 The lattice in this case is just a theoretical framework for defining the regularized 
model, and not a very convenient numerical tool. 
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However, the free energy density of U(N) is a gauge invariant quantity 
and its expansion in 1/N should not depend on the particular gauge. For 
the purpose of 1/A-counting, background gauges [22] in fact are much easier 
to use than covariant MAG. Contrary to MAG one cannot define the BRST- 
symmetry of background gauges on the lattice [2H] since the lattice gauge 
group is compact j20] • But these are renormalizable gauges that are well 
defined to all orders in perturbation theoryj22J 121] ■ This suffices for our 
purpose. Background gauges and MAG share the crucial properties that the 
background ga^ does not renormalize |22t "2""fj and that it couples minimally 
to the fields. Since background gauges are linear, they do not constrain the 
color flow and do not change the 1/N counting of a diagram. 

The background gauge in our case is defined by a gauge-fixing part of the 
Lagrangian of the form, 

= ^ \P*yl*\ - C b a V,{V,C)t . (19) 

T>n and in Eq.(|19jl are, respectively, the background covariant deriva- 
tive (with the connection ga^ defined in Eq.(j2}) and the ordinary covariant 
derivative (with connection gV^). C and C denote the ghost and anti-ghost 
fields and a is the gauge parameter. Upon shifting the gauge field by the 
constant and Abelian background a M , the premise that all time derivatives 
occur as background covariant derivatives holds in these gauges. 

Apart from rigorously defining perturbative propagators and introducing 
a set of adjoint ghost fields, there are no constraints on the color summations 
in the background gauge fixing of Eq.([T""]). These gauges therefore do not 
modify any of the previous arguments with regard to the order in N of a 
perturbative diagram. 

4.2 Regularization and Renormalization 

Background gauges are renormalizable to all orders in perturbation theory [2*1] . 
We nevertheless have to show that the previous scaling argument is not 
spoiled by the renormalization procedure. Although the free energy den- 
sity superficially diverges quartically, the specific heat at constant volume 
(Cy = —Td^F) is only quadratically divergent. Cy may, for instance, be 
regularized by analytic continuation in spatial dimensions only. 

The free energy density is recovered by integration of the specific heat 
with the boundary conditions that the specific entropy, —8tF, and the free 
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energy density, F, vanish at T = 0. This is equivalent to subtracting from 
the free energy density any contribution that is linear in the temperature. 
For D = 3 — £ spatial dimensions, a dimensionally regularized perturbative 
contribution to the specific heat is of the form, 

C V (T, N, g 2 - e, fi) = Td 2 T G(T, N, g 2 - e, //) , (20) 

where G(T, N, g 2 ; e) is the formal expression of the vacuum graph in D spatial 
dimensions, g 2 = g 2 [i~ e = X/N is the renormalized dimensionless coupling 
and /i is the renormalization scale. The diagrammatic argument of section 3.1 
implies that the contributions of planar gluonic vacuum graphs in the center- 
symmetric background depend on T and N in the particular combination, 

G S2 (T,N,g 2 ;pL,e) = NTf S2 (T/N,\;e,fM) = ^fs 2 (k Nfi/T, e) . (21) 

The subtraction of a constant term and of a term proportional to T from 
G amounts to the subtraction from / of a term proportional to N/T and 
of a T- independent constant. Possibly divergent terms from planar vacuum 
diagrams that are proportional to N 2 and N thus do not contribute to the 
specific heat nor to the free energy density. 

Further, since the free energy is a physical quantity, /(A; iV/x/T, e) does 
not depend on the renormalization point /i. In the renormalization scheme 
(RS), /(A; Nfi/T, e — > + ) therefore is a function of the renormalization 
group invariant effective coupling X(T/A RS ) only. 

The free energies of center-symmetric planar U(N) for different iV are 
proportional only if the temperature is measured in terms of an asymptotic 
scale parameter, Ars, that does not depend on N. To determine this finite 
renormalization, it is sufficient to, for instance, demand that the deconfine- 
ment temperature Td(N) of planar U(N) be the same for all N. The scaled 
free energy N 2 F S2 (T, N) then does not depend on iV at any temperature 
below T d . 

Eq. (j21j) shows that the coupling X(Njj) is a function of Nfi rather than 
of the renormalization point fi only. Large values of iV correspond to large 
values of fi - and to weak coupling. For large N this implies that, 

24tt 2 

X(iiN) . (22) 

Eq.p^ suggests that the confining phase can be explored perturbatively 
at sufficiently large N. This weak-coupling confinement regime was first 
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noticed by Polchinski[ll] while exploring the analogy between string theory 
and large-N gauge theory. However, the background of Eq.(JU) is expected 
to be unstable for temperatures T > T^. Setting the renormalization mass 
fi ~ Td in Eq. (J22j) . the unstable regime corresponds to couplings A < A^ 
with 

UN) ~ 24 fl x • (23) 

For any finite value of N, the phase transition occurs at a (perhaps small) 
but nevertheless finite value of the coupling. An asymptotic perturbative 
expansion thus is not possible in the center-symmetric phase for any fixed 
value of N. Due to Eq.(J23J), a perturbative evaluation of the (leading) O(N ) 
contribution to the free energy could nevertheless be reasonably accurate. 

Since the usual 1/iV-expansion of SU(N) is algebraic in l/N, the loga- 
rithmic dependence of the coupling on N in Eq.(|23|) is somewhat unexpected. 
However, the center-symmetric orbit of Eq.fjHJ) is described by a connection 
that is itself iV-dependent. This leads to a non-trivial iV-dependence of the 
momentum scale in planar diagrams. The usual ultraviolet behavior of the 
model then results in a logarithmic dependence on iV of the effective cou- 
pling. The perturbative analysis of U (N) gauge theory in the confining phase 
in this sense becomes self-consistent at large N. 

5 Other contributions to the free energy of 
center-symmetric U(N) gauge theory 

We saw that the contribution from planar 't Hooft diagrams to the free energy 
in the center-symmetric phase is of order l/N 2 only. 'T Hooft diagrams with 
the topology of a torus may superficially contribute to the free energy density 
in order N°. To conclude that the free energy of U(N) indeed is of order N° 
in a center-symmetric l/N expansion we have to consider some remaining 
possibilities. 

5.1 No contributions to the free energy of O(N) 

Fields in the fundamental representation of the group explicitly break the 
center symmetry and superficially could give rise to contributions to the free 
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energy that are of order N. We will see that there in fact are no such contri- 
butions in an expansion about the center-symmetric background of Eq.(J3J. 
The argument is rather similar to the one employed in the gluonic case. Vac- 
uum diagrams that superficially are of order N are planar diagrams with 
one fundamental color loop only [A sphere with a hole, topologically a disc 
D 2 .] A typical 't Hooft diagram of this kind is shown in fig. 2b. We now 
have that L p = L c and can augment to composite loop momenta as before. 
The difference to planar gluonic vacuum diagrams is the absence of an extra 
perimeter loop over color only. This suppresses such contributions by a fac- 
tor of iV compared to the planar gluonic ones of fig. 2a. The sums over the 
time-components of composite loop momenta now extend over half-integer 
multiples of the fundamental frequency 2nT/N. The previous scaling argu- 
ment shows that such vacuum diagrams contribute to the free energy density 
in order l/A 3 : 

F D2 (T,g 2 ,N) = Tfn 2 (T/N,X) = ^f D2 (\(T/A)) (24) 

Below we explicitly find that this is also true for contributions of order A . 
There thus are no contributions of order N 2 or N in the expansion of the 
free energy density of a U(N) gauge theory about the center-symmetric back- 
ground of Eq.(j4*]). Since planar contributions to the free energy from adjoint 
and fundamental fields vanish in the limit of large N, the center-symmetric 
planar U (N) model approaches a topological theory without dynamical de- 
grees of freedom. 



5.2 The free energy of center-symmetric U(N) gauge 
theory to order A 

The previous diagrammatic analysis does not extend to the 1-loop contribu- 
tion to the free energy density. We explicitly compute it for an U (N) gauge 
theory with Np Dirac fermions in the fundamental representation. The rel- 
evant quadratic part of the Lagrangian is, 

N r 1 1 

6 

v a) 



o,6=l ' 
N F N 

+(V»C a b )(V,C b a )] + E [*£rA*i + ™M*j] ■ (25) 
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In Eq. (j25|) the 7 M are the hermitian Euclidean Dirac matrices that satisfy 
Ifilu + lul^i = 1<^ii/ • The time component of the background covariant deriva- 
tive T>i for the fundamental and adjoint representation is given in Eq.© 
and Eq.(J7J) respectively (V^ = for spatial indices /i ^ 4). The gluon- (V^) 
and ghost- (C, C) fields satisfy periodic boundary conditions in temporal 
direction whereas the fermions (\F,-, tyj) are anti-periodic. 

Since the free energy does not depend on the gauge parameter, we may 
for simplicity choose the Feynman-like gauge a = 1 to compute it. For the 
constant background of Eq.(JU), the eigenvalues of the operator V^V^ are 
readily obtained and the functional integral over quadratic fluctuations about 
this background can be formally performed. For D = 3—e spatial dimensions, 
the regulated contribution to d^F (T, N; /i, e), of the non-interacting model 
is: 



d*F (T,N;n,e) 



2^a=l J (2tt)" 



(26) 



EL 2 In 



k' 2 +(2nT/N) 2 (nN+a-b) 2 
/,2 



-4 Eft ^ 



k 2 + {2TrT/N) 2 (nN+a-l/2) 2 +m 2 



Noting that nN + a ranges over all integers, Eq.()26|) simplifies to, 



k 2 + {2nT/N) 2 n 2 
,,2 



k 2 +(2irT/N) 2 (n-l/2) 2 +m 2 



(27) 



This expression converges for D < 1 spatial dimensions and thus is at most 
quadratically divergent. Scale invariance of the non-interacting model de- 
fined by Eq. (|23j) implies the absence of quadratic divergences in the massless 
casej2S|. 

The contribution to the free energy of a non-interacting massless bosonic 
degree of freedom at temperature T is finite and for D = 3 spatial dimensions 
is[26j, 

T 4 7T 2 

F boS on(T, m = 0) = — — . (28) 

yu 

That from a non-interacting massive fermionic degree of freedom is finite as 
well 



f ermion 



(T,m) 



;„ 2 7' 2 g (-1 



2tt 2 



n=l 



rr 



-Kzinm/T) . 



(29) 
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In Eq. (l29j) K 2 (z) is the K-Bessel function normalized so that for small argu- 
ments K 2 (\z\ ~ 0) = 2/z 2 . [Note that |F boson (T, m = 0) = F fermion (T,m = 
0) < Ff erm i on (T,m) < 0. The last inequality results because z 2 K 2 {z) is a 
monotonically decreasing function of its argument on the positive real axis, 
with z 2 K 2 (z) < 2 for all z > 0. The contribution of massive fermions to the 
free energy density is exponentially small for T <C m.} 

With the integration conditions that the free energy density and the spe- 
cific entropy vanish at zero temperature, ( F (0,N) = <9t-F o (0,A0 — 0), the 
specific heat completely specifies the free energy density. One can read off the 
1-loop contribution to the free energy density from Eq.()27p: to lowest order in 
the coupling, the free energy density of center-symmetric U (N) gauge theory 
at temperature T is that of 2N 2 non-interacting bosonic degrees of freedom 
and 4NNp (massive) fermionic degrees of freedom but at a temperature of 
T/N. Using Eq.flEJ and Eq.(J2H one has, 

F (T, N) = 2N 2 F boson (T/N, m = 0) + AN J2?=i F fermton (T/N, m 3 ) 

= -IS + 3m Ef=i J2n=i ^ (^) 2 K 2 {nNm 3 /T) . (30) 

We thus find the same behavior in iV for the 1-loop contribution to the 
free energy density as for planar diagrams. We in fact used the same ar- 
guments, combining color- and momentum- flow to a single composite mo- 
mentum / 4 (n, a) = (2nT/N)(nN + a). Although the center-symmetric back- 
ground effectively leads to an ^-dependent rescaling of the temperature, it 
is perhaps more appealing to view this as an (for sufficiently large N) almost 
complete cancellation of individual contributions to the free energy (see ap- 
pendix A). Since background of Eq.(jlJ) essentially shifts frequencies by a 
fraction of the fundamental frequency 2irT a partial cancellation can occur: 
for non-interacting fields this phase shift effectively amounts to a change in 
boundary conditions and the free energy density is sensitive to this change 5 . 
However, it is remarkable that this cancellation is almost complete in planar 
U(N) at large N. 

5 The contribution to the free energy density of a bosonic degree of freedom satisfying 
anti-periodic boundary conditions (corresponding to a shift of the frequency by ttT) for 
instance is positive. 
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5.3 Contributions to the free energy of order N° 

Center-symmetric planar U(N) in the limit iV — » oo is devoid of physical 
degrees of freedom and a topological model. Fortunately, the scaling argu- 
ments we used to show the absence of contributions to the free energy of 
orders N 2 and N break down when there are more independent momentum- 
than color- loops, that is when L p > L c . One then has at least one loop 
momentum that cannot be augmented to a composite momentum that in- 
cludes the color flow. Such contributions to the free energy density do not 
scale with N. The corresponding 't Hooft diagrams either include more than 
one fermion loop or are non-planar. The leading contribution of O(N ) is 
given by diagrams that topologically either are a torus or a sphere with two 
holes [i.e. a disk with a hole]. These two classes of 't Hooft diagrams corre- 
spond to contributions to the free energy of order N° from non-interacting, 
colorless, asymptotic glueball- and meson-states respectively. These are the 
stable asymptotic states^] at large N. 

The fact that the free energy is O(N ) for sufficiently large N implies 
that in the center-symmetric background of Eq.pjl all asymptotic states are 
color singlets. No higher dimensional representation of the global color group 
contributes. The center-symmetric expansion "confines" color in this sense. 
[Note that this definition of confinement is somewhat stronger than Wilson's 
screening criterion for static color charges in the fundamental representation. 
The latter cannot be applied in the presence of light meson states.] 

To conclude that U(N) gauge theory confines color charge at sufficiently 
large N one would have to show that the background is stable against fluctua- 
tions in some (low) temperature regime. The fact that the strong coupling ex- 
pansion of lattice gauge theory confines and is center-symmetric suggests that 
this might be the case at sufficiently large effective coupling A(T/A) > X d . To 
further conclude that the more realistic S77(3)-model confines color charge 
at low temperatures one in addition has to show that there is no (deconfin- 
ing) phase transition at some finite N > 3. Neither of these issues will be 
discussed further here. Let us instead look at some interesting aspects of the 
previous analysis. 
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6 Infrared-Finite Perturbation Theory at Fi- 
nite Temperature 

The suppression of planar contributions to the free energy density shows that 
color momentum is essential. For a background like Eq.(j5J) that breaks the 
center symmetry (and corresponds to vanishing color momentum £), Linde 
observed [T] that the perturbation series of a non-Abelian gauge theory is 
infrared divergent at any finite temperature. The most infrared divergent 
vacuum diagrams are all planar and superficially are of order N 2 . The center- 
symmetric background of Eq. (J3J) provides an effective infrared cutoff of order 
2itT/N for all the coset excitations. Since it appears via the covariant deriva- 
tive, this infrared cutoff is not entirely equivalent to an effective gluon mass. 
Unlike an effective mass, it does not regulate the Abelian sector of the model 
in the infrared. The gauge bosons of an Abelian ?7(l) Ar -model on the other 
hand do not interact directly, and the infrared behavior of such models is 
regular when all charged fields are massive [27]. Even though the option 
of massive off-diagonal fields is not available for a U(N) gauge theory, the 
center-symmetric background of Eq. (JIJ provides an infrared cutoff that works 
rather similarly: it shifts the infrared singularity of any coset field propagator 
from k 2 = to k 2 = -4ir 2 T 2 j 2 / N 2 for some integer N/2 > j > 0. Note that 
although some of these "masses" are rather small at large N, they do not 
depend on the coupling A. 

In the center-symmetric background of Eq.pjl the perturbation series 
thus is free of infrared divergences without resummation. Of course, when 
the effective coupling is sufficiently small (at high temperatures) this back- 
ground presumably is not stable (see below). The center-symmetric back- 
ground of Eq.(J3} cures the pervasive infrared problem of the perturbative 
expansion at low temperatures only. Although the perturbation series may 
not converge in this regime, its mere existence to all orders does define the 
model formally. The regularization of perturbative infrared divergences by 
the center-symmetric background, however, reshuffles contributions to the 
1/N expansion of the free energy. It does so in a manner that is consistent 
with confinement in this phase. 
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7 Stability and (Veneziano's) Vector Ghosts 



The result that planar U(N) gauge theory practically has no degrees of free- 
dom at large N, implies that center-symmetric planar SU(N), although de- 
void of colored asymptotic states, is not a thermodynamically stable model. 
Center-symmetric SU(N) gauge theory nevertheless can be a perfectly good 
physical model because the subset of planar diagrams does not give the lead- 
ing contribution in 1/N. There then is no reason why this subset of diagrams 
should define a thermodynamically viable physical model. Planar diagrams 
are generated by the Cuntz algebra .29] rather than by a bosonic or fermionic 
one. There is no proof that such a field theory is thermodynamically stable. 

The instability of center-symmetric planar SU(N) follows immediately 
from the previous result for the planar U(N) model without fundamental 
fields. The color singlet "photon" decouples in this case and the free energy 
density of U(N) is just that of the corresponding SU (N) model and of a free 
photon. Eq. (j28|) together with the previous result for U(N) implies that the 
free energy density of center-symmetric planar SU(N) is, 

C W (T)~ + 0(l/iV 2 ). (31) 

The positive contribution to the free energy of center-symmetric SU(N) is 
of O(N ) and can be interpreted as due to a massless, color-singlet vector 
ghost that compensates the degrees of freedom of the massless, color-singlet 
"photon" of center-symmetric planar U(N). 

Veneziano |30j has shown that a massless color-singlet vector ghost in pla- 
nar gluonic SU(N) could saturate the axial Ward Identities and solve the 
L^(l)-problem at large N. Eq. (J31j) is evidence for the existence of a vector- 
ghost in the confining phase of the planar model. Whether this vector-ghost 
couples to the axial current in the manner Veneziano suggests, cannot be de- 
termined from the free energy. To have a viable confining phase, the vector- 
ghost of the planar model would have to either decouple by itself (as all states 
in the planar truncation do) or be part of a BRST-quartet j2HI that does not 
contribute to the free energy. 

As discussed in the previous sections, there are additional contributions 
to the free energy density of O(N ) in the center-symmetric phase that are 
described by non-planar diagrams. These non-planar contributions to the 
free energy density depend on the effective coupling A. It is at least con- 
ceivable that massless bound states form when A > A^ that complete the 
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BRST quartet and compensate the contribution to the free energy of the 
vector ghost. Since the vanishing of ghost contributions to the free energy is 
necessary for the stability of a center-symmetric phase, the critical coupling 
at which this occurs is a lower bound for A^. 

The fact that the non-interacting (Casimir) part of the free energy density 
of SU(N) is positive (from Eq.flHJ) with Eq.fll|) for all N >2, implies that 
the center-symmetric phase is not stable at small effective coupling A(T/A) ~ 
0. This is consistent with the expectation that the center symmetry is broken 
for T > T d > 0. 

When non-planar contributions to the free energy of O(N ) are included, 
SU(N) could be thermodynamically stable in the center-symmetric phase at 
sufficiently large effective coupling A(T/A) > A^. 

8 Discussion and Conclusion 

The configuration of Eq.(0J) is an absolute minimum of the classical action 
of unbroken U(N) gauge theory. It is invariant under the discrete global 
Zn center-symmetry of the Yang-Mills action. The description of a center- 
symmetric orbit of vanishing curvature is unique in the infinite volume limit 
of a periodic lattice. It is unique even at finite volume if the spatial topology 
is that of a three-sphere, since the only non-contractible Wilson loops in this 
case are the Polyakov loops in temporal direction. 

Although the background of Eq. (0J) is an absolute minimum of the clas- 
sical Yang-Mills action, it is not an absolute minimum of the free energy 
density at all temperatures. Due to its symmetries, this orbit is always an 
extremum of the free energy, but to lowest order of perturbation theory this 
extremum is a maximum. To the extent that higher order perturbative cor- 
rections are negligible, the center symmetry is broken for sufficiently small 
effective coupling, that is at sufficiently high temperatures [TH]- The strong- 
coupling expansion of lattice gauge theory suggests that a center-symmetric 
phase is thermodynamically preferred at low temperatures when the effective 
coupling is sufficiently strong. 

An expansion in 1/N could be an appealing alternative to the strong 
coupling expansion in this non-perturbative regime. The center-symmetry of 
the gluonic sector is preserved in every order of the perturbative expansion 
about the background configuration of Eq.(jlJ). Certain qualitative conclu- 
sions about the 1/A^-expansion of the free energy can be obtained by exam- 
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ining this perturbative series. There are no contributions to the free energy 
of order N 2 or N at any perturbative order. For large N, a center-symmetric 
planar truncation of U (N) approaches a topological field theory without dy- 
namical degrees of freedom. This confirms the result of Gocksch and Neri[S] 
that the free energy of lattice gauge theory in the planar limit does not de- 
pend on the temperature (and therefore vanishes) at large N. The leading 
contribution is of order iV° and due to vacuum diagrams that could represent 
the free energy of color singlet quark-antiquark mesons and glueballs. 

This is as one expects for the confining phase of a U(N) gauge theory. 
Perhaps more significant is that perturbative calculations at finite tempera- 
ture in principle are feasible in the background of Eq. (@J . The severe infrared 
divergences of ordinary perturbation theory observed byjTj do not occur in 
this center-symmetric expansion. The reduced coupling A furthermore be- 
comes weak in the confining phase for sufficiently large N (see section 4.2). 
Some aspects of confinement therefore may be accessible in a perturbative 
framework The perturbative analysis of the model in the vicinity of 
the center-symmetric minimal action orbit in this sense is (self) consistent. 
However, the existence of a phase transition at a perhaps very small but nev- 
ertheless finite value of the coupling restricts the accuracy of a perturbative 
analysis in the confining phase at any fixed value of N. More quantitative 
results perhaps can be expected from summing classes of perturbative con- 
tributions. 

Planar SU(N) turns out to be thermodynamically unstable at sufficiently 
large N due to a massless color-singlet vector-ghost. This is not an artefact 
of our treatment but rather a consequence of the fact that the [/(l)-photon 
decouples and the free energy of planar U (N) in the confining phase vanishes 
as iV — ► oo. We speculate that the vector ghost couples to the axial current 
in the manner conjectured by Veneziano |3Uj. It then should be part of a 
BRST-quartet |2H1 that, as a whole, does not contribute to the free energy. 
This could be the case at temperatures T < Ta when other contributions to 
the free energy density of O(N ) are included. 

Although the center symmetry of the action is explicitly broken by fields 
in the fundamental representation, they do not contribute in order N to 
the free energy in an expansion about the center-symmetric background 
of Eq.pjl. Spontaneous chiral symmetry breaking could perhaps be inves- 
tigated in this background: to lowest order in the coupling and for large N 
several fermionic degrees of freedom are almost zero modes in the background 
of Eq.pjl. Whether this is sufficient to trigger spontaneous chiral symmetry 
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breaking has not been explored. If so, the thermodynamic instability of the 
background for T > would imply that the chiral- and deconfinement- 
transition temperatures coincide. 
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A Casimir contributions to the free energy of 
center-symmetric U(N) 

We here calculate contributions to the free energy density in the center- 
symmetric background of Eq.(j3J) to zeroth order in perturbation theory with- 
out recourse to the scaling argument. The calculation explicitly shows that 
contributions of individual degrees of freedom cancel. 

Consider the regularized expression 6 for the free energy density given 
by Eq.([26|). The free energy density to lowest order in the coupling can be 
decomposed, 

N N F N 

F (T, N) = 2 J ( T > (a " h )l N i °) " 4 E E J ( T > ( a " l l 2 )l N \ m ^ ' ( 32 ) 

a,b=l j=\ a=l 

into individual contributions I(T, 5; m) that depend on the phase 5 and mass 
m associated with a particular degree of freedom. J(T, 5; m) is formally given 
by, 

I(T,5;m)= lim £ V) / ln(fc 2 + m 2 + 4vr 2 T 2 (n + 6) 2 ) (33) 



6 For simplicity and to easily include fermions, the following computation uses dimen- 
sional regularization. The gluonic contribution to the free energy of SU (N) has also been 
computed in lattice reerularizationpTT] . Evaluating Neuberger's lattice result confirms that 
this gluonic contribution to the free energy is 0(1/ N 2 ) at the center-symmetric background 
and sub-leading for sufficiently large N. 
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The expression of Eq. (J33|) is not well defined. However, making use of the 
fact that the free energy density F and the specific entropy —dxF vanish at 
T = 0, it suffices to obtain the second derivative d\l. Explicitly taking the 
derivatives in Eq.([26|h the expression for <9fJ(T, 8; m) may be written, 

«2r^x \ ^ fd D k^~ D 4n 2 T(n + 8) 2 

^ /(T ' 5;m) = zfe- 2. J {2 n)° k 2 + m 2 + 4n 2 T 2 (n + 8) 2 (34) 

n=— oo v ' y ' 

2(k 2 + m 2 ) 



x 



1 + 



k 2 + m 2 + An 2 T 2 (n + 8) 2 _ 

Km f) 2 T ^ °^ f dDk ^ D ^-A(fc 2 +m 2 +47r 2 T 2 (n+5) 2 ) 



[Note that the final (finite) result does not depend on the renormalization 
point fM. The latter was introduced to have a free energy density with the 
canonical dimension.] We next evaluate the momentum integrals in the last 
expression of Eq.(|34|). 



To separate the summation over the integers n from the dependency on 
8 it is convenient to use the Fourier-representation of the Gaussian : 



Tipn 



E e -47r 2 AT 2 (n+5) 2 = f°° dp c - p 2/^T 2 X) c 2TTi P 8 

n 

-n 2 /(4T 2 A) (36) 




— L= ^1 + 2 y cos(27™5)e-" 2 /( 4T2A ) 



[Note that the second expression for the sum is dual to the original one in the 
sense that the "radius" of the temporal direction has been inverted 4AT 2 — > 
1 / \4ir 2 XT 2 ) .] Substituting the last expression for the sum into Eq.([35|) and 
noting that the constant term in the braces of Eq. (j36|) does not survive 
differentiation, one finds that in D = 3 dimensions, 



CO 



I(T, 8;m) = — > cos(27rn5) 



^ g -Am 2 -n 2 /(4T 2 A) 
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E 



71=1 



cos(27rn<5) / Tm 



2tt 2 



n 



K2(nm/T) . 



(37) 



The integration constants in Eq.(|37J) have been determined so that the free 
energy density and the specific entropy vanish at T = 0. 

The result of Eq.(|37J) can be checked in various limits: the free energy 
density of a non-interacting bosonic degree of freedom satisfying periodic 
boundary conditions is obtained with 5 = 0; 5 = 1/2 corresponds to the 
(positive) free energy of a non-interacting bosonic degree of freedom satisfy- 
ing anti-periodic boundary conditions, etc. One evidently can achieve some 
cancellation in the total free energy by mixing several degrees of freedom sat- 
isfying different boundary conditions. In the center-symmetric background 
of Eq.(jlJ), the gluons effectively satisfy different boundary conditions. If the 
fundamental and anti-fundamental color indices of the gluon (in the back- 
ground gauge) are a and b, the corresponding shift of the Matsubara fre- 
quency is 27rT(a — b)/N. Note that diagonal Abelian degrees of freedom 
(with a = b) do not suffer a phase shift and that all bosonic degrees of 
freedom are iV-periodic. 

The gluonic part of the free energy density in Eq. (|3*2*|) is found by evalu- 
ating, 



N 



a,b=l 



N 



2 cos(27m(a - b)/N) = 



,2nin(a-b)/N 



a,b=l 
N 



a=l 



t 2irian/N 



+ C.C. 



(38) 



Here 5 n ^N is the Kronecker symbol (one when n is an integer multiple of N 
and zero otherwise). Semiclassically, contributions to the free energy density 
from gluonic paths with (n mod N) ^ windings thus cancel completely. 
Using Eq.(|3*7|) (for m — >• 0) and Eq. (|3*5|) the gluonic contribution in Eq.| 



is. 



N 



N 



2^J(T,(a-p;0) = -2T 4 £ £ 



cos(2vrn(a - b)/N) 



a,b=l 



=1 n=l 



7T% 4 
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riSr 2N 2 T 4 tt 2 

n 2 (kN) 4 ~ 45A 2 ' ( ' 

This verifies the scaling argument for the first term of Eq. ()30)) . The fermionic 
contribution in Eq. (j30)l is similarily obtained from Eq. (j37j) and Eq. (j32j) . We 

have that, 

N 1 N 

^cos(27m(a- 1/2) /N) = ^ ^ e 2 ™ (a " 1/2)/iV + c.c. 

a=l a=l 

= Nj2(-V k Sn,kN ■ (40) 

k 

Using Eq. (jlUj) and Eq. (j3*7j) in Eq. (|3^j) gives the second term of Eq. (|3T?|) . 
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